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基于 矢量 相似 性 的 多 元 滤波 方法 研究 
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摘 要 : 为 了 避免 彩色 图 像 滤波 时 排序 难题 ， 在 探究 失 量 相似 性 的 基础 上 ， 提 出 一 种 针对 彩色 图 像 的 多 元 滤波 方法 。 
首先 ， 在 RGB 彩色 空间 中 ， 利 用 彩色 和 拓 量 间距 离 和 角度 定义 其 相似 性 测度 ， 以 刻画 与 人 类 视觉 感知 相 一 致 的 彩色 相 
似 性 ; 其 次 ， 以 上 述 彩色 相似 性 为 准则 设计 并 构建 一 种 彩色 多 元 滤波 方法 ， 并 就 其 相关 参数 对 滤波 性 能 的 影响 进行 了 
深入 分 析 和 研究 ; 最 后 ， 为 了 验证 提出 方法 的 有 效 性 ， 将 其 应 用 于 标准 彩色 影像 和 彩色 逐 感 影像 滤波 。 实 验 结果 可 知 ， 
提出 的 方法 不 仅 有 效 解 决 了 传统 滤波 方法 的 排序 难题 ， 而 且 克 服 了 因 滤 波 使 图 像 模糊 ， 边 缘 不 清 等 问题 ; 另外 ， 将 提 
出 的 多 元 滤波 与 传统 方法 进行 实验 对 比 , 结果 表明 其 不 仅 能 有 效 滤 除 多 种 类 型 噪声 , 而 且 较 好 地 保持 了 原始 图 像 信息 ， 
使 图 像 信息 保 真 清晰 ， 其 视觉 效果 优 于 传统 ， 客 观 评价 指标 也 有 较 大 幅度 的 改善 和 提升 。 
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(OO Research on multivariate filtering method based on vector similarity 
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Abstract: In order to avoid the problem of sorting in color image filtering, this paper defines and proposes a multivariate filter 
method for color image based on the study of vector similarity. Firstly, in RGB color space, the similarity measure between the 


color vectors ls defined by the distance and the angle to describe the color vector similarity consistent with the human visual 


和 perception. Secondly, using the color similarity criterion, the color multivariate filter is designed and constructed, and the 
influence of the related parameters on the filtering performance is analyzed and studied in depth. Finally, in order to verify the 
effectiveness of the proposed method, it is applied to the actual remote sensing image filtering. It can be seen from the 
experiment that the proposed method not only effectively solves the problem of sorting the traditional filtering methods, but 
also overcomes the problems of blurring the image and unclear edges due to the filtering. In addition, compared with the 
results of other traditional methods, the results show that it not only can effectively filter all kinds of noise, but also maintains 
the original image information better, making the image information clear and fidelity. The visual effect of this filtering 
method is superior to the traditional, its objective evaluation index also has greatly improved. 
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现 


方法 具有 很 强 的 针对 性 ， 图 像 和 噪声 的 特性 不 同 ， 采 用 的 滤波 
方法 亦 有 所 不 同 。 从 图 像 表 达 形 式 来 看 ， 图 像 滤波 可 分 为 灰 度 
图 像 滤波 和 彩色 图 像 滤 波 两 大 类 。 灰 度 图 像 滤波 以 灰 度 值 作为 
图 像 质量 及 其 处 理 结果 精度 贴 。 由 于 产生 原因 不 同 ， 常 见 的 唯一 依据 ， 实 现 较 为 容易 。 与 灰 度 图 像 不 同 ， 彩 色 图 像 是 多 元 
像 噪声 包括 高 斯 、 椒 盐 等 多 种 类 型 。 图 像 滤波 则 是 以 图 像 保 真 像 , 涉及 多 个 颜色 通道 , 而 且 各 通道 间 存 在 着 强烈 的 相关 性 。 
为 目标 的 图 像 处 理 技术 ， 甚 作用 在 于 从 被 噪声 污染 的 图 像 中 恢 “ ”因此 ， 灰 度 图 像 滤 波 方 法 不 能 直接 移植 于 彩色 图 像 ， 也 使 彩色 


噪声 产生 于 图 像 采 集 、 传 输 和 处 理 的 各 个 环节 ， 直 接 影 


于 


| 


| 


复原 始 图 像 信息 ， 同 时 保持 图 像 中 目标 边缘 和 细节 信息 ， 因 此 ” 图像 滤波 的 实现 更 复杂 、 更 困难 。 
在 图 像 处 理 中 有 着 十 分 重要 的 实际 意义 。 一 般 而 言 ， 图 像 滤 波 对 于 彩色 图 像 ， 噪 声 不 仅 使 其 颜色 失真 、 清 晰 度 下 降 ， 而 
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实际 


上 且 也 直接 影响 图 像 的 处 理 及 应 用 所 ， 针 对 噪声 常用 的 滤波 方法 
有 标量 和 矢量 滤波 法 。 标 量 滤 波 法 是 对 彩色 图 像 的 三 个 颜色 通 
道 分 别 进行 滤波 ,然后 合成 各 通道 滤波 结果 得 到 彩色 滤波 图 像 。 
这 类 方法 没有 考虑 彩色 图 像 中 三 个 颜色 分 量 间 的 有 机 联系 ， 可 
能 产生 原 图 像 所 没有 的 颜色 ， 进 一 步 造 成 图 像 失真 ， 甚 至 出 现 
局 部 不 协调 的 色彩 ， 从 而 破坏 了 色调 和 边缘 细节 信息 。 矢 量 滤 
波 法 是 将 彩色 像素 作为 彩色 空间 矢量 来 处 理 ， 充 分 考虑 各 颜色 


常 可 


] 时 ， 矢 量 排序 标准 不 同 ， 则 滤波 输出 也 有 
成 YMF、VDF、DDF 等 非 线性 滤波 方法 09-24。 
E 离 作为 排序 依据 时 ， 则 窗口 中 每 


所 有 矢量 距离 的 累加 和 4; 可 表示 为 


此 规则 获得 排序 序列 为 


通道 分 量 之 间 的 联系 ， 在 消除 噪声 、 保 持 色调 、 保 护 边 缘 与 细 
节 上 有 较 好 的 鲁 棒 性 。 因 此 ， 彩 色 图 像 的 矢量 滤波 法 较 标 量 滤 
波 法 更 合理 、 更 有 效 ， 其中， 矢量 中 值 滤 波 (vector median filter, 
VMF)、 矢 量 方向 滤波 (vector direction filter VDF)、 矢 量 方向 距 
离 滤波 (vector direction distance filter，DDF) 等 方法 受到 普遍 关 
注 S 引 ,但 这 些 方法 在 滤波 时 均 采用 最 大 数量 平滑 , 经 常会 因 损 
坏 图 像 的 边缘 和 细节 信息 而 导致 图 像 失 真 。 除 了 上 述 矢量 滤波 
方法 外 ，Lukac 等 人 外 提出 选择 加 权 的 矢量 中 值 滤波 ; 武 昆 等 
人 60 将 四 元 数理 论 引 入 到 彩色 图 像 滤波 领域 ， 孙 建 召 等 人 吕 
利用 模糊 神经 网 络 进行 图 像 噪声 滤波 等 。 这 些 改 进 的 滤波 方法 
使 其 性 能 在 某 些 方面 得 到 一 定 的 提升 和 改善 ， 但 由 于 彩色 图 像 
表示 形式 的 特殊 性 及 排序 等 问题 ， 所 以 均 存 在 一 定 的 局 限 性 2 


为 了 避免 彩色 图 像 滤波 时 排序 的 难题 ， 本 文 以 彩色 空间 矢 
量 相似 性 为 准则 ， 设 计 并 构建 了 一 种 针对 彩色 图 像 的 多 元 滤波 
方法 。 该 方法 不 仅 消除 了 因 排 序 带 来 的 滤波 问题 ， 而 且 能 较 好 
滤 除 图 像 噪声 ， 并 保持 原始 图 像 细 节 特 征 ， 也 不 会 因 产 生 新 的 
颜色 而 造成 图 像 失真 。 与 传统 的 滤波 方法 比 对 ， 其 性 能 得 到 了 
较 大 的 改善 和 提升 。 


1 “理论 基础 


给 定 彩色 图 像 z={zi(xi, y), Oo eeP, ie7 } 为 定义 在 图 像 域 
P 上 的 彩色 影像 ， 其中, T= {1, 2, … n}，(xi, 为 像素 i 的 平面 
坐标 , i 为 像素 索引 , n 为 像素 个 数 , P 为 图 像 域 , zi= (zin, zic, zi8) 
为 像素 i 在 RGB 中 所 对 应 的 彩色 矢量 ，zz、zic、zi 分 别 表示 
彩色 矢量 zi 的 红 、 绿 、 蓝 分 量 。 以 像素 i 为 中 心 的 滤波 窗口 彩 
色 矢量 集合 Wi= (zi, zi) = {zj, jeJ}， 其 中 ，ieJ，i'eNi，Ni 为 像 
素 i 的 邻 域 像素 集合 ，j 为 滤波 窗口 中 像素 索引 ,J = {1, 2, .…， 
m}，m 为 滤波 窗口 中 矢量 数目 。 彩 色 图 像 矢 量 滤波 就 是 按 一 定 
的 规则 和 顺序 逐个 获取 滤波 窗口 并 实现 滤波 输出 ， 这 样 彩色 图 
像 滤波 就 转化 成 彩色 空间 矢量 的 滤波 问题 04481。 
1.1 典型 矢量 滤波 
在 进行 矢量 滤波 时 ， 将 以 像素 ;为 中 心 的 滤波 窗口 W: 中 以 
个 彩色 矢量 ， 按 照 一 定 的 规则 进行 排序 并 构成 一 个 升序 序列 : 


T 


Zz 


2 Ly SS jeJ 


(1) 
其 中 , zw 表示 第 j 矢量 的 排序 统计 量 ， 则 排序 最 小 矢量 za 
和 最 大 矢量 zm 就 分 别 是 该 滤波 窗口 的 最 小 和 最 大 滤波 输出 。 


去 除 在 矢量 距离 


zj 到 所 有 矢 


则 式 (3) 所 对 应 式 (了 ) 序 列 


4,= Dd(z,,z) jke 
k=1 


A(1)<AQ)...< AO)...<Am) 


(2) 


) 是 当前 滤波 窗 


距离 累加 和 最 小 自 
值 滤波 VMF,， 这 种 滤波 方法 输 


其 作为 滤波 输出 ， 即 就 是 天 


出 为 窗口 内 原 有 像素 , 能 


4 噪声 像素 。 


4 和 角度 来 进行 排序 时 ， 则 每 个 矢量 


则 式 (5) 对 应 式 (1) 序列 


三 


的 累加 和 OG; 可 表示 为 
O),= > ,0,2) jkel 
= 


类 似 的 排序 序列 为 
OW <O0) ...<O0 ...<Ou) 
Pp 的 zw 是 当前 滤波 窗口 


矢量 
素 。 


外 考虑 矢量 相关 性 的 特征 ， 


vo 
汽 
ee 
区 


有 时 会 将 矢量 


(4) 


矢 


作为 滤波 输出 ， 即 就 是 所 谓 的 
它 能 有 效 地 去 除 色调 相差 较 大 的 噪声 像 


距离 和 相对 方向 同时 作为 排序 依据 ， 则 滤波 窗口 中 每 个 矢量 和 


到 所 有 矢量 的 距离 与 角度 累加 和 的 乘积 Bj 表示 为 


也 发 4 


以 距离 和 相对 方向 


这 时 ， 式 (7) 所 对 应 式 (1) 序 列 
所 有 矢量 距离 与 角度 累加 和 乘积 最 小 的 矢量 ， 这 就 是 距离 方向 


B,=> 0,z) dz,z) jkel 
k=1 天 =1 


度 〉 作 为 矢量 排序 依据 
BO <B 2) ...<B 0) ...<B un) 


的 zw 是 当前 滤波 窗 


(6) 


对 ， 可 得 到 


(7) 


到 


滤波 DDF， 既 能 去 
周 相差 较 大 的 噪声 像素 。 


除 矢量 距离 相差 较 大 的 


心 


色 i 


ES 


虹 声 像素 ， 又 能 消除 


均 为 窗口 内 原 有 像素 ， 不 会 引入 新 


的 


E 了 变化 ， 从 而 致使 彩 1 


息 在 滤波 过 程 中 受到 损伤 和 影 


ay 


存在 
算 量 随 滤波 窗口 包 


平滑 ， 没 有 被 噪声 污染 的 像素 


际 应 用 


通用 
像 、 保 留 
1.2 矢 
役 而 言 ， 彩 色 图 像 中 未 受 噪声 影响 的 同 质 


有 较 高 的 相似 性 ， 而 受到 噪声 影响 后 其 相似 性 就 会 降低 ， 本 文 


曾 大 呈 几 何 速度 增长 ， 导 致 难 
， 在 实际 应 用 时 


方面 均 存 在 一 定 的 局 限 性 。 


,图 像 中 的 轮 廊 、 边 缘 、 细 节 等 
， 在 抑制 噪声 和 1 方 本 
如 时 均 要 进行 大 量 矢量 排序 ， 运 
以 满足 一 些 实 

| ， 其 实 难以 找到 一 种 合适 、 

以 这 些 方法 在 去 除 还 原 图 
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正 是 利用 
的 相似 性 程度 ， 六 
效 避 免 和 解决 了 彩 


这 个 特征 ， 


以 此 作为 彩色 图 像 多 元 滤波 的 基础 ， 这 也 
色 矢 量 排序 的 难题 。 


引入 相似 性 测度 的 概念 以 刻画 彩色 矢量 间 


相似 性 测度 是 研究 和 分 析 目 标 间 关系 的 一 种 十 分 重要 的 量 
度 方法 ， 在 图 像 处 理 、 信 息 检 索 、 模 式 识别 、 数 据 挖 气 、 多 元 
数据 分 析 、 图 像 检索 等 中 都 有 着 十 分 广泛 的 应 用 。 

对 于 任 一 像素 i 所 构成 的 滤波 窗口 彩色 矢量 集合 Wi = {zj 
jeJ}， 给 定 彩 色 矢量 对 z, zre Wi，j,keJ， 则 其 矢量 相似 性 测度 
可 定义 为 

HZ) =e "cos (KO(z,,21)) (8) 


其 中 ，kie[0,%) 和 ke[0,1] 为 控制 相似 程度 的 参数 ，d(zj, zt) 和 
Q(zj, zp e[0, x/2) 分 别 为 两 彩色 矢量 间 的 距离 和 夹 角 ， 可 表示 为 


3 Up 
(Seal| je 
g=1 


de 9) 
max||zn -zu zl lz -zi |] p=% 
> zc 
cosgGcjzo)= 王 0 (10) 
(Br) (ar) 
其 中 : p 为 矢量 范式 ， 在 这 里 取 p=2， 即 dz z 为 矢量 间 的 欧 
几 里 德 距离 。 由 式 (8) 可 以 看 出 ， 定 义 的 矢量 相似 性 测度 取 值 区 


何 晓 军 ， 等 


多 元 滤波 的 实现 奠定 了 基础 。 


4 
当 d=1,0=m4 时 kl、k2 对 相似 性 测度 的 影响 


图 1 

然 ， 式 (8) 中 参数 h 决定 了 两 个 矢量 空间 距离 上 的 相似 
性 ， 当 参数 hh = 0 时， 矢量 相似 性 仅 与 角度 有 关 ; 同 理 ， 参 数 
忆 决 定 了 两 个 矢量 在 角度 上 的 相似 程度 ， 当 有 b= 0 时 ， 矢 量 相 
似 性 仅 与 距离 有 关 。 


I 


2 ”算法 描述 

彩色 图 像 z 的 多 元 滤波 就 是 按照 顺序 对 其 每 一 像素 点 i 所 
对 应 的 滤波 窗口 Wi 进行 处 理 ， 并 实现 滤波 矢量 输出 vy:， 即 yi;= 
£(Wi)， 当 前 像素 i 在 窗口 中 对 应 的 矢量 设 定 为 zz?。 具 体 实 现 过 
程 如 下 。 


间 为 [0, 1]。 这 样 ， 彩 色 图 像 中 两 个 像素 点 颜色 在 视觉 感知 上 越 据 式 (8)-(10) 计 算 矢 量 间 的 相似 性 测度 并 形成 一 个 WixWi 
相近 时 ， 其 对 应 到 RGB 空间 的 彩色 矢量 在 距离 和 角度 上 就 越 ”相似 性 关系 矩阵 Ri: 
接近 ， 即 其 距离 和 夹 角 越 小 ， 相 似 性 测度 的 值 越 大 ， 则 表示 两 Ri= [UL zd)] j,key, iel (11) 
者 具有 较 高 相似 性 ; 反之 ， 当 两 个 像素 点 颜色 在 视觉 感知 上 差 由 式 (11) 得 到 的 相似 性 关系 矩阵 Ri 满足 以 下 两 个 条 件 : 
别 较 大 时 , 其 矢量 间距 离 或 夹 角 较 大 , 其 相似 性 测度 的 值 越 小 ， a) 自 反 性 : yw,z)=1,， jeJ 
则 表示 相似 程度 较 低 。 当 然 ， 当 是 同一 颜色 时 ， 两 个 彩色 矢量 b) 对 称 性 : Li zb =4(zw 2),， j,keJ 
重合 ( 即 为 同一 矢量 )， 距 离 及 夹 角 均 为 0， 相 似 性 测度 为 1， 则 对 于 zie Wi 的 平均 相似 性 测度 可 表示 为 
则 表征 其 相似 程度 最 高 。 彩 色 图 像 的 实际 视觉 感知 均 与 矢量 相 ] 之 . 
=— PAT) jkel (12) 
似 性 相符 合 。 ma 
一 般 而 言 ， 人 类 对 彩色 认 知 存在 差异 性 ， 相 似 性 会 因 人 的 式 (12) 中 平均 相似 性 测度 4; 表示 矢量 zj 与 Wi 中 所 有 矢量 
视觉 感知 而 不 同 ， 另 外 ， 彩 色 图 像 的 边缘 、 细 节 信 息 等 均 具有 的 相似 性 测度 的 平均 值 ， 由 这 些 平均 相似 性 测度 中 最 大 值 所 对 
定 的 非 确定 性 ， 使 得 人 们 很 难 精 确 地 区 分 这 些 边缘 细节 像素 应 的 矢量 构成 一 个 矢量 集合 Ui 可 表示 为 
和 噪声 , 所 以 式 (89) 中 引入 参数 后 和 局 ,通过 改变 其 值 能 达到 人 = {zreWi, 4= max{Aj}} jtey,iel (13) 
为 调整 两 个 矢量 的 相似 性 程度 的 目的 ， 即 对 于 同样 两 个 矢量 ， 由 式 (13) 即 可 获取 Wi 中 所 有 矢量 相似 性 测度 平均 值 最 大 的 
于 参数 hh、 忆 选择 的 不 同 ， 则 具有 不 同 的 相似 性 测度 值 。 如 ”矢量 ， 也 就 是 当前 窗口 中 最 “中 心 ” 矢 量 。 这 样 ， 即 可 得 到 一 
车 dlzj,z) =1、9(zj,z) =z/4， 两 个 矢量 间 相似 性 测度 变化 ， 如 。 ”个 与 最 “中 心 ”矢量 对 应 的 相似 性 聚 类 及 对 应 的 相似 性 测度 ; 
图 1 所 示 , 从 图 中 能 够 清楚 地 看 出 , 针对 同样 的 两 个 彩色 矢量 ， ”根据 实际 滤波 要 求 ， 如 果 给 定 一 个 阔 值 a(0 <a<1)， 则 由 此 相 
当 参 数 后 和 e 取 不 同 值 时 ， 由 式 (8) 所 得 到 的 模糊 相似 性 测度 。” 似 性 聚 类 中 矢量 相似 性 测度 大 于 给 定 闵 值 a 的 矢量 构成 一 个 割 
不 同 ， 这 与 人 类 视觉 感知 的 模糊 性 和 不 确定 性 相 一 致 。 通 过 实 。 集 Wsc WW， 称 为 a 市 集 ， 即 
验 分 析 可 知 , 参数 后 和 妨 的 不 同 组 合 直接 决定 了 两 个 矢量 的 相 Wia= {zp u(t) > ao} tey (14) 
似 性 程度 ， 随 着 h、 厂 的 增加 ， 矢 量 间 的 相似 性 测度 减 小 ， 相 则 在 当前 滤波 窗口 W 中, 基于 矢量 相似 性 的 多 元 滤波 输出 
似 性 降低 ; 反之 ,相似 性 升 高 ; 参数 后 比 握 对 相似 性 测度 的 影 ”可 表示 为 
响 更 大 ， 相 似 性 测度 对 于 后 的 变化 更 加 敏感 。 上 述 相似 性 测度 z, zeW 
义 能 很 好 的 刻画 两 个 矢量 间 的 相似 性 关系 ， 这 也 为 彩色 图 像 "二 z EW 3 
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从 式 (15) 可 知 ， 
即 原 为 原 当前 
当前 滤波 窗口 中 最 
“中 心 ”矢量 


即 形成 滤波 后 彩 


f『 矢量 zl; 


均 为 当前 
十 于 给 定 彩色 
色 图 像 了 = 


当前 矢量 在 Wis 中 时 ， 其 滤波 窗口 输出 vi 
当前 矢量 zi 不 在 Wis 中 时 ， 其 输出 yi 为 
“中 心 ”矢量 zw， 而 无 论 是 当前 矢量 还 是 最 
滤波 窗口 中 的 原 图 像 像 素 矢 量 。 
图 像 z= {zi, ie 1 了}， 通 过 上 


述 操 作 #， 
{vi ie7T}， 其 中 , 1={1,2,..., n}。 


这 种 处 理 不 仅 能 够 更 好 的 保留 原 图 像 的 细节 特征 ， 同 时 也 避免 


因 矢 量 排序 带 


和 提升 。 


3 ”实验 结果 与 分 析 
3.1 多 元 滤波 窗口 中 相关 量 的 获取 


来 的 滤波 问题 ， 算 法 效率 及 性 能 有 了 较 大 的 改善 


滤波 窗口 是 多 元 滤波 的 基本 单元 和 区 域 ， 为 了 进一步 分 析 


其 中 各 个 量 的 获取 和 确定 。 假 设 彩色 图 像 z 中 ， 


以 某 一 像素 i 


为 中 心 构造 一 个 3x3 像素 的 滤波 窗口 Wi, 其 对 应 的 矢量 集 Wi= 


{zp jed}s d= 42 0 


9}， 九 个 彩色 矢量 分 别 为 : 


z1=(35, 47, 


49)、 z2=(232, 236, 236)、 z3=(85, 97, 99)、 z4=(29, 34, 13)、 zs=(143, 


表 1 滤波 窗口 中 矢量 相似 性 测度 


和 zo=(12, 18, 41)， 


和 矩阵 Ri 及 
1 可 以 得 出 ， 平 均 相 
的 矢量 集合 为 U;= {zs}， 


窗口 中 最 “中 心 
发 现 一 个 最 优 矢量 作为 输出 


何 晓 军 ， 等 


Chingiy 会 作 期 刊 


目 似 1 羽 久 合集 波 方法 研究 


145, 147) 、z6=(9, 23, 45)、2z7=(143, 137, 146) 、zs=(56, 59, 70) 


其 中 ， 当 前 滤 


波 像素 i 对 应 的 矢量 为 z1。 


为 了 简化 实验 复杂 度 ， 选 取 参 数 有 = 0.001， 厂 = 0.2， 利 用 
式 (8)~(12) 计 算 矢量 间 的 相似 性 测度 并 形成 一 个 WixWi 相似 性 
F 均 相似 性 测度 Aj, 计算 结果 如 表 1 所 示 。 由 表 


似 性 测度 4 的 最 大 值 为 0.9124， 其 所 对 应 


即 zs 是 当前 窗口 最 “中 心 ”的 矢量 ， 


前 滤波 矢量 zi 未 包含 在 w- 割 集中 ， 
”的 矢量 zs， 这 就 说 明 此 方法 在 允许 范围 内 能 
Hh， 而 其 仍 为 原 图 像 矢 量 。 


此 矢量 所 对 应 相似 性 聚 类 的 相似 性 测度 分 别 为 :{0.9683, 0.7409， 
0.9455, 0.9553, 0.8653, 0.9343, 0.8698, 1.000, 0.9326} 。 当 滤波 阔 
值 a=0.9 时 ， 其 w- 割 集 Wia= { zu za,z4 z6 zs z9}， 因 当前 滤波 
矢量 zi 在 a- 制 集中 ， 所 以 滤波 窗口 输出 vi; 即 为 原 矢量 xz1， 这 充 
分 体现 了 此 滤波 方法 在 允许 范围 内 尽量 保留 原 窗口 及 像素 的 原 
有 特征 ;， 当 滤波 阀 值 a= 0.97 时 ， 其 w- 割 集 Wia= {zs}， 此 时 当 


所 以 滤波 输出 产 即 为 此 滤波 


矢量 2Z2 Z3 Z4 ZX6 Z7 Z8 Z9 

RGB (35,47,49) (232,236,236) (85,97,99) (29,34,13) (143,145,147) (9,23,45) (143,137,146) (56,59,70) (12,18,41) 
Zl 1.0000 0.7180 0.9169 0.9845 0.8386 0.9618 0.8426 0.9683 0.9599 0.9101 
Z2 0.7180 1.0000 0.7832 0.7076 0.8561 0.6920 0.8515 0.7409 0.6903 0.7822 
Z3 0.9169 0.7832 1.0000 0.9035 0.9145 0.8839 0.9187 0.9455 0.8816 0.9053 
Z4 0.9845 0.7076 0.9035 1.0000 0.8265 0.9748 0.8307 0.9553 0.9748 0.9064 
Z5 0.8386 0.8561 0.9145 0.8265 1.0000 0.8080 0.9920 0.8653 0.8061 0.8786 
Z6 0.9618 0.6920 0.8839 0.9748 0.8080 1.0000 0.8121 0.9343 0.9927 0.8955 
Z7 0.8426 0.8515 0.9187 0.8307 0.9920 0.8121 1.0000 0.8698 0.8104 0.8809 
Z8 0.9683 0.7409 0.9455 0.9553 0.8653 0.9343 0.8698 1.0000 0.9326 0.9124 
Z9 0.9599 0.6903 0.8816 0.9748 0.8061 0.9927 0.9326 1.0000 0.8943 


0.8104 


3.2 ”多 元 滤波 实验 


为 了 直观 订 


150x150 像素 的 经 典 彩色 图 


E 价 这 种 多 元 滤波 方法 的 性 能 ， 分 别 选 取 
像 Pepper"” 和 分 辨 率 为 1 m, 150x150 


像素 的 Ikonos 高 分 辨 率 遥 感 影像 作为 测试 图 像 ， 如 图 2(al)(a2) 


所 示 ; 分 别 加 入 强度 为 0.08 的 椒盐 噪声 构成 噪声 图 像 ， 如 图 
2(b1)(b2) 所 示 ; 式 (8) 中 计算 矢量 相似 性 测度 时 ， 仍 然 选取 参数 
太 = 0.001, 有 a= 0.2, 滤波 窗 


0.97。 经 过 多 元 滤波 后 的 输 


以 看 出 ， 此 多 元 滤波 方法 
的 去 除 噪声 , 而 且 很 大 程度 
未 因 滤 波 造成 图 


像 信 


为 3x3 像素 的 矩形 区 域 , 闵 值 a = 


出 图 像 ， 如 图 2(c1)(c2) 所 示 ， 从 原 
图 像 (图 2(al)(a2)) 与 滤波 后 图 像 (图 2(c1)(c2)) 直 观 对 比 效 果 可 


应 用 奠定 了 基础 。 


x 


有 很 好 的 滤波 效果 ， 不 但 能 够 较 好 


保留 了 实际 影像 的 边缘 及 细节 信息 ， 
息 模糊 


， 避 人 免 失 真 ， 为 彩色 图 像 的 进一步 


图 2 图 像 滤波 
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3.3 ”参数 对 滤波 性 能 影响 的 分 析 与 评价 通过 前 面 的 算法 描述 可 以 看 出 ， 因 多 元 滤波 方法 以 矢量 相 

实际 评价 滤波 方法 效果 及 性 能 主要 是 通过 主观 和 客观 两 个 ” 似 性 为 基础 ， 所 以 式 (8) 中 引入 的 参数 后 、 心 及 式 (14) 中 的 阔 值 
方面 。 主 观 评价 是 通过 人 的 主观 感觉 对 图 像 质量 的 优 劣 做 出 评 ”a 直接 影响 滤波 性 能 。 其 中 参数 hh、ke 决定 了 两 个 矢量 之 间 的 
定 ， 上 述 实 验 即 通过 此 方式 进行 评价 ， 该 方式 具有 很 大 的 主观 。 相似 程度 , 当 k=0 时 , 两 个 矢量 间 的 相似 程度 与 距离 没有 关系 ， 
性 ， 受 不 同 观察 者 、 不 同 图 像 类 型 和 环境 的 影响 ， 使 其 具有 一 。” 仪 与 相对 角度 有 关 ， 即 滤波 就 是 以 矢量 角度 作为 依据 进行 ， 即 
定 的 不 确定 性 。 因 此 ， 在 实际 图 像 滤 波 质量 评价 中 主要 以 客观 。 变 为 VDF 滤波 ; 同 理 ， 当 k=0 时 ， 就 变 为 VMF 滤波 ， 即 两 个 
评价 为 主 ， 客 观 评价 就 是 通过 量化 指标 进行 滤波 图 像 的 评价 。 矢量 间 的 相似 性 与 相对 角度 没有 关系 ， 仅 与 距离 有 关 ， 仅 以 距 
般 而 言 , 图 像 质量 特征 不 同 , 客观 评价 的 标准 也 会 有 所 不 同 。 ”高 作为 滤波 依据 ; 闵 值 a 决定 了 滤波 窗口 输出 当前 像素 点 的 概 
但 在 实际 应 用 中 ， 由 于 很 难 评估 滤波 方法 消除 噪声 以 及 对 未 受 率 。 
噪声 影响 像素 的 保持 能 力 ， 所 以 对 其 性 能 的 评价 是 一 件 十 分 复 首先 , 评价 参数 有 在 不 同 噪声 干扰 下 对 滤波 效果 及 性 能 的 
杂 的 工作 。 本 文采 用 应 用 最 为 广泛 的 标准 均 方 差 normalized 影响。 假设 已 = 0.2， 阔 值 wc= 0.5， 而 后 从 0 变化 到 10 时 ， 多 
mean square error NMSE)， 通 过 统计 两 个 数字 图 像 之 间 像 素 值 ”元 滤波 的 性 能 NMSE 评价 曲线 ， 如 图 4(a)-(d) 所 示 。 从 图 4 可 
分 布 的 不 同和 差异 来 进行 多 元 滤波 方法 的 性 能 评价 。 若 zi 表示 。” 知 ， 当 如 = 0.001 时 ， 多 元 滤波 方法 对 于 噪声 模式 1 的 滤波 效 
原 彩色 图 像 中 像素 i 所 对 应 的 矢量 ， 而 六 表示 其 估计 值 ， 即 滤 。 果 最 好 ， 如 图 4(a) 所 示 ; 当 可 = 0.02 时 ， 其 对 于 噪声 模式 2-4 


波 后 该 像素 的 具体 值 ， 则 NMSE 可 表示 为 的 滤波 效果 最 好 ， 如 图 4(b) 所 示 。 从 参数 五 整体 影响 的 变化 趋 
Sz vl 势 看 ， 多 元 滤波 针对 噪声 模式 3 椒盐 噪声 的 滤波 性 能 尤为 
NMSE = 一 一 (16) 显著 。 
2 za or 
分 析 参 数 对 滤波 性 能 的 影响 实验 ， 选 取 150x150 像素 的 四 
“Pepper" 作 为 测试 图 像 ， 如 图 2(aD) 所 示 ， 在 这 里 分 别 用 高 斯 噪 。 和 ” 


声 、 椒 盐 噪 声 及 两 者 结合 进行 污染 干 捧 ， 噪 声 类 型 及 相关 参数 
如 表 2 所 示 。 


表 2 噪声 类 型 
模式 噪声 类 型 
1 高 斯 (方差 0.005) i 
2 高 斯 (方差 0.005) + 椒盐 (强度 0.02) 有 
3 椒盐 (强度 0.02) 
4 椒盐 (强度 0.02) + 高 斯 (方差 0.005) 
原 测试 图 像 经 过 上 表 品 声 进行 干扰 后 的 噪声 图 像 ， 如 图 3 
所 示 ， 图 3(a)~(d) 分 别 对 应 表 2 中 四 种 噪声 模式 产生 的 噪声 图 图 4 参数 kl 对 噪声 模式 1-4 滤波 性 能 的 NMSE 评价 


同 理 , 评价 参数 忆 对 滤波 性 能 影响 时 , 设 定 参数 操 = 0.001， 
闷 值 a = 0.5， 参 数 be[0,1]，NMSE 变化 曲线 ， 如 图 5 所 示 ， 
从 图 示 变 化 可 以 看 出 ， 当 ha = 0.2 时 对 所 有 噪声 模式 其 滤波 效 
果 及 性 能 达到 最 优 。 


-er Model 1 
一 本 Model2 
一 ~- Model 3 
-EH- Model14 


NMSE 


(Gd 
图 3 噪声 图 像 (a) 模式 1 (b) 模 式 2 (c) 模 式 3 (d) 模 式 4 图 5 参数 对 噪声 模式 1-4 滤波 性 能 的 NMSE 评价 
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一 本 Miodel 2 
一 | Model 3 
-EModel14 


NMSE 


6 浆 值 a 对 噪声 模式 1-4 滤波 性 能 的 NMSE 评价 

另外 , 闵 值 a 是 用 于 调节 滤波 时 保留 原 像 素 的 概率 ， 当 有 
= 0.001， 有 如 = 0.02， 闪 值 a 在 0~1 变化 时 ，NMSE 曲线 如 图 6 
所 示 。 图 中 变化 情况 表明 ， 当 a<0.6 时 ， 对 于 所 有 噪声 模式 ， 
滤波 效果 最 好 。 
基于 上 述 实验 及 分 析 ， 可 以 看 出 ， 参 数 后 、 刀 及 阔 值 x 均 
直接 影响 多 元 滤波 的 性 能 ， 所 以 在 实际 中 一 般 要 根据 具体 情况 
进行 确定 。 
3.4 性 能 对 比 实 验 
下 面 将 本 文 多 元 滤波 方法 与 矢量 中 值 滤波 VMEF、 矢 量 方向 
滤波 VDF、 距 离 方向 滤波 DDF 等 典型 滤波 方法 进行 对 比 实验 ， 
进一步 研究 和 分 析 其 性 能 。 根 据 前 面 的 实验 结果 分 析 ， 本 文 多 
元 滤波 方法 参数 后 = 0.001、 忆 = 0.2、a = 0.5, 滤波 窗口 均 为 3x3 
像素 ， 四 种 测试 噪声 图 像 如 图 3 所 示 。 

针对 噪声 模式 1 的 VMF、VDF、DDF 及 本 文 方法 滤波 结 
果 分 别 如 图 7(al)~(a4) 所 示 ; 针对 噪声 模式 2、3、4 的 滤波 结 
果 如 图 7(b1-b4)~(d1-d4) 所 示 。 从 滤波 效果 对 比 可 以 看 ， 本 文 方 
法 (图 7(a4)~(d4)) 有 具有 良好 的 滤波 效果 及 性 能 , 能 较 好 的 保留 原 
图 像 的 细节 特征 ， 且 能 使 图 像 边缘 光滑 。 

当然 ， 四 种 方法 的 滤波 性 能 及 效果 有 时 通过 主观 视觉 很 难 
准确 进行 评判 , 在 这 里 , 仍然 采用 上 述 标准 均 方 差 (NMSE) 来 进 
行 评价 和 分 析 ， 其 滤波 性 能 NMSE 曲线 如 图 8 所 示 ， 其 中 图 
8(a) 为 四 种 滤波 方法 针对 噪声 模式 1 的 滤波 性 能 NMSE 评价 ， 
即 对 应 图 7(al)-(a4); 同 理 ， 图 8(b)-(d) 为 针对 噪声 模式 2、3、4 
的 NMSE 评价 曲线 ， 且 分 别 对 应 图 7(bl-b4)-(d1-d4)。 从 NMSE 
线 可 以 清楚 的 看 出 ， 本 文 提出 的 多 元 滤波 方法 与 传统 的 
VMF、VDF、DDF 等 方法 相 比 ， 其 客观 评价 指标 有 较 大 的 改善 
和 提高 ， 针 对 上 述 四 种 噪声 均 具有 较 好 的 滤波 效果 及 性 能 。 


4 ”结束 语 


为 了 解决 彩色 图 像 滤波 问题 ， 文 章 在 RGB 彩色 空间 ， 从 
矢量 相似 性 角度 出 发 ， 利 用 距离 和 角度 共同 定义 相似 性 测度 ， 
以 刻画 和 描述 矢量 间 与 人 类 视觉 感知 相 一 致 的 相似 性 关系 ; 
以 此 为 基础 和 准则 ， 构 造 并 提出 了 一 种 针对 彩色 图 像 的 多 元 滤 
波 方法 。 该 方法 不 仅 能 避免 滤波 时 矢量 排序 的 难题 ， 而 且 具 有 
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保留 原始 像素 和 去 除 噪声 像素 2 
像 边缘 及 细节 特性 。 通 过 实验 可 知 其 具有 较 好 的 滤波 效果 ， 尤 
其 是 当 参 数 语 = 0.001， 妨 = 0.2， 阔 值 a<0.6 时 ， 该 方法 滤波 性 
能 较 优 。 另 外 ， 文 章 还 将 其 与 常用 滤波 方法 进行 实验 对 比 ， 通 
过 客观 评价 指标 NMSE 的 变化 曲线 表明 , 该 方法 能 有 效 的 去 除 
多 种 类 型 的 噪声 ， 能 较 好 地 保持 原 图 像 信息 ， 其 视觉 效果 及 评 
价 指标 均 优 于 传统 滤波 方法 ， 具 有 一 定 的 推广 和 实用 价值 。 


间 的 折衷 能 力 ， 较 好 保持 了 图 


图 7 针对 四 种 类 型 噪声 模式 的 滤波 结果 


Wm 


本 (d) 


图 8 不 同 滤波 方法 针对 不 同 噪声 模式 的 NMSE 评价 
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